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ABSTRACT 



Aims. Implement a matched filter (MF) cross-correlation algorithm in multipole space and compare it to the standard Angular Cross 
Power Spectrum (ACPS) method. Apply both methods on a Integrated Sachs Wolfe (ISW) - Large Scale Structure (LSS) cross corre- 
lation scenario and study how sky masks influence the multipole range where the cross correlation signal arises and its comparison to 
theoretical predictions. 

Methods. The MF requires the inversion of a multipole covariance matrix that, under non-full sky coverage < 1), is gener- 
ally non-diagonal and singular. We choose a Singular Value Decomposition (SVD) approach that enables the identification of those 
modes carrying most of the information from those more likely to introduce numerical noise, (that are dropped from the analysis). 
We compare the MF to the ACPS in ISW-LSS Monte Carlo simulations, focusing on the effect that a limited sky coverage has on the 
cross-correlation results. 

Results. Within the data model s = t + am where t is Gaussian noise and m is a known filter, we find that the MF performs compara- 
tively better than the ACPS for smaller values of and scale dependent (non-Poissonian) noise fields. In the context of ISW studies 
both methods are comparable, although the MF performs slightly more sensitively under more restrictive masks (smaller values of 
fsky). A preliminary analytical study of the ISW-LSS cross correlation signal to noise (S/N) ratio shows that most of it should be 
found in the very large scales (50% of the S/N at / < 10, 90% at / < 40 - 50), and this is confirmed by Monte Carlo simulations. 
The statistical significance of our cross-correlation statistics reaches its maximum when considering / 6 [2, l„,cx], with /„,„,- 6 [5, 40] 
for all values of observed, despite of the smoothing and power aliasing that aggressive masks introduce in Fourier space. This 
/-confinement of the ISW-LSS cross correlation should enable a safe distinction from other secondary effects arising at smaller (higher 
/-s) angular scales. 

Key words. (Cosmology) : cosmic microwave background. Large Scale Structure of the Universe 



1. Introduction 

Auto and cross-correlation analyses are crucial in the study 
of the Cosmic Microwave Background (CMB) anisotropics. 
This is due not only to the fact that the theory can only 
predict statistical properties of the intrinsic intensity and 
polarization anisotropics (and hence auto-correlation tests must 
be conducted in order to compare theory to observations, see 
|Hu & Dodelson(2002)| for a review), but also due to the presence 
of secondary anisotropics and foreground emission that add 
up to the measurements in the microwave range. These other 
components must be identified and separated from the intrinsic 
ones generated at the surface of last scattering, and therefore 
cross-correlation analyses to other data sets probing the sources 
of this secondary emission must be carried out. This has been 
done practically for all CMB experiments, from COBE data 
dSmoot et al.( 1992)1 [Bennett et s£(199^ all the way to WMAP 
data ( [Bennett et al.(2003)| [Spergel et al.(2007)| l. These cross 
correlation techniques may be either based in real space (like the 
angular two point correlation function), in Fourier space (like 
the Auto and Cross Angular Power Spectrum), or in wavelet 
space ( [Cayon et al.(2000)| [Larson et al.(2005)| . 



space the statistical covariance matrices between diff'erent modes 
are particularly simple, and so is the comparison of theory to ob- 
servations. It is in this space where theoretical expectations for 
other secondary effects present in the CMB are also displayed, 
and where the constraints on the cosmological parameters are set 
(e.g., [Dunkley et al.(2008)| [Reichardtet al.(2008)| l. However, 
there are two practical issues that tend to complicate this theory 
to data comparison: the presence of Cosmic Variance in the 
large angular scales (that is, the sample variance due to having 
only one single sky to look at) and the coupling of different 
Fourier modes whenever not the entire sky is subject to analysis 
(as it happens in practice for current and future CMB and LSS 
surveys like ACT (Kosowsky(2003)| l, SPT ( [Ruhl et al.(2004)l l, 
DUN^ SNAlQ etc). These two eff'ects are of particular 
relevance in the study of the Integrated Sachs-Wolfe (ISW) 
eff'ect ( [Crittenden & Turok (1996)| l: the ISW arises in the large 
angular scales, and since its frequency dependence is identical 
to that of the intrinsic CMB fluctuations, it must be identified 
via cross-correlation tests to Large Scale Structure surveys that 
are likely to cover only a fraction of the sky. 

In this work we generalize the matched filter cross- 
correlation method to multipole space in the context of CMB 



In the linear theory that characterizes the intensity and po- 
larization anisotropics of the CMB, predictions are done in the 
Fourier space of the 2D sphere, that is, in multipole space. In this 
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Studies. We compare it to the standard Angular Cross Power 
Spectrum in different scenarios, and show that the former is ei- 
ther equivalent or superior to the latter We also perform this 
comparison in Monte Carlo simulations of the ISW effect, with 
similar results. The method is developed in Section (2), whereas 
a first comparison to the Angular Cross Power Spectrum is given 
in Section (3). A detailed analysis of the signal to noise ratio 
of ISW cross-correlation measurements is provided in Section 
(4), where the matched filter method is again compared to the 
Angular Cross Power Spectrum. Finally, in Section (5) we dis- 
cuss our results and conclude. 



2. The matched filter method 

2.1. A brief description 

Our first goal is to estimate the level of presence of some known 
signal m in some measured data array s, which is therefore de- 
composed as s = t -t- am. We shall assume that t is a Gaussian 
vector (which will be regarded as noise) whose covariance ma- 
trix C is known. Given the Gaussian assumption, C completely 
characterizes t. As shown in, e.g., |G6rski et al.(I996)[ the mini- 
mization of the quantity 



= - am)iiC ')/,/s - 0-111)^, 



yields the following estimates for a and its formal error: 



t'C Im 
m'C im' 



1 



m'C im 



(1) 



(2) 



Note that the superscript f denotes transpose. The difficulty 
usually lies in the inversion of the covariance matrix for 
long data arrays t and/or for close-to-singular covariance 
matrices C. The first scenario was already addressed in 



Hernandez-Monteagudo et al.(2006) where this technique was 
appUed in separated subsets of data, and then the covariance 
among different subsets was computed separately. Here, we 
shall also consider the case where C is singular or close to 
singular. 

Indeed, the use of the matched filter is very ex- 
tended in CMB analyses (e.g., |Rubino-Martm et al.(2000)"l 



Hernandez-Monteagudo & Rubifio-Martm(2004) 



Hernandez-Monteagudo et al.(2004)| [Hansen et al.(2005)| ), 
but it has been mostly restricted to real space. In works like 
that of [Hansen et al.(2005)| it was also implemented in Fourier 
(multipole) space, but only after approximating the covariance 
matrix as diagonal, assumption that we shall avoid here. 



2.2. The covariance matrix in multipole space 

We will focus our analyses on real signals defined on 2D spheres. 
These are usually decomposed on an spherical harmonic basis as 
follows: 



s(fi) = ^ ^ a;,,„y/,„,(n), 

/=/,„„, m=-l 



(3) 



with n denoting a direction on the sky (or a position in the 
sphere). If s is real, then the multipole coefficients verify a; - 
(-l)'"a*i^, with the symbol "*" denoting complex conjugate. 
This limits the number of degrees of freedom per / to 2/ -i- 1. 



The m - multipole is by definition real, so the 21 remain- 
ing degrees of freedom can be assigned to the real and imagi- 
nary parts of the o; ,,, multipoles with magnetic number ranging 
from m - 1 to m = /. I.e., for a given multipole /, the multipole 
array fl/_,„ will be decomposed into all + I dimensioned array 
(uL ui , uL vL .., vi), where uL = ai o, u\ , u', contain the real 



parts of ai,„=i i, and Vj, v{ the imaginary ones. Since we will 



simultaneously consider all multipoles / e [/„„„, /„ 
the multipole array 

a = (uo, u, v). 



x], we define 



(4) 



where 

Uo = (m|;"' 

u = (m'™' 
and 

v = (y'™",. 



'„„■„+ 1 



(5) 
(6) 

(7) 



The dimension of u and v is given by ni 2 - ImaxQmax + 
l)/2 + l„,ax -i- 1 - {Iminilmm " 1) + hnin), whcrcas the dimen- 
sion of Uo is simply «;j - l„,ax - Imin + 1, SO the total dimen- 
sion of a reads «/ = n/j -H 2 X n^- If s(n) is an isotropic, 
Gaussian distributed signal over the whole sphere (fsity - 1), 
then the correlation matrix of the coefficients is diagonal: 
{aimiav ,m')*) - Ci6ij'6m,m'- (Notc that due to isotropy there is 
no dependence on m). Likewise, we have that (C)ij = {aiaj) 
is diagonal in such case. This fact makes the inversion of the 
covariance matrix in equation ([Tj) trivial. Let us now relax the 
assumption on having s(n) defined over the full sphere. In an as- 
trophysical context, if some parts of the sky are lacking data, i.e., 
if fsk\ < 1, the covariance matrix C will no longer be diagonal 
and for large enough l„,ax it will also be singular A traditional 
matrix inversion is likely to either fail or provide inaccurate re- 
sults. (Note that the accuracy of the inversion can be tested by 
running Monte Carlo simulations and comparing the dispersion 
of the recovered a's with the actual prediction of equation (|2])). 

In these circumstances, we perform a S VD decomposition of 
the covariance matrix. 



C = R'AR, 



(8) 



where A is a diagonal matrix (containing the eigenvalues of C) 
and R is a rotation orthogonal matrix (R'R = I). Note that since 
C is symmetric and positive definite, the eigenvalues should all 
be positiv^l The SVD decomposition sorts the eigenvectors 
according to the magnitude of the eigenvalues, so the first 
eigenvectors are those containing more information about C, 
whereas the last ones are the most likely to introduce numerical 
noise. Note that there is always some numerical error in our 
estimates of the covariance matrix, since it is computed through 
a finite number of Monte Carlo realizations (10,000 in this 
work). Therefore, this decomposition provides a way to rotate 
the vector a into its principal modes, and permits distinguishing 
those having most of the information from those being more 
affected by numerical error (which can be safely projected 
out of the analysis). In practice, we neglected all eigenvectors 
whose eigenvalues were smaller that a given fraction e of the 
first (largest) eigenvalue. (For most cases, the choice e = 10"^ 
yielded optimal results). The inversion of C after the SVD 
decomposition becomes straightforward, enabling an easy 



^ In practice, we find that for dense and close to singular cases of C, 
some eigenvalues (of small absolute value) were negative. 
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implementation of the matched filter as given by equation (|2]i. 
Note that, unlike in |Gorski( 1 994)| or [Mortlock et al.(2002)| we 
are not worried in building a new set of orthonormal functions 
in the patch of the sky under analysis, nor we attempt to perform 
component separation ( [Bouchet et al.(1999)| l. In all those works, 
the techniques used were in some way close to ours, but their 
goals were different. 

As we shall see below, we may be interested in applying the 
matched filter in different /-bins. One can readily find that, given 
the outcome of the matched filter in two different Z-bins Op and 
aq, their covariance is given by 



m'r T r 'm 

"V 99 ? 

{m!pC-plmp){m'qCqqmq) ' 



(9) 



Here, Cpp and Cqq denote the covariance matrices of the noise 
t for Z-bins p and q, respectively, whereas Cpq is the covariance 
matrix for the noise in different I binfl: Cpq - (tptq). For a set 
of /-bins we shall obtain a vector of measured a^-s, whose com 
bined;^f^ will be given by 



Y^&piC ^)p,q&q 



(10) 



An overall detection level for a given set of /-bins and as will 
be provided by this statistic. Another statistic providing the 
level of detection is the variance weighted average for the lip's, 
here defined as y6: 



(11) 



We will show below that, in ISW studies, the distribution of 
the Up-s will be very close to Gaussian, and therefore Gaussian 
will also be the distribution of p. (Note that the matched filter 
method, as defined from a minimization of the statistic given in 
equation ([TJ, is only optimal if the noise is actually Gaussian 
distributed). The diagonal elements of the matrix C can be com- 
puted via equation (|9]l or via numerical simulations: the agree- 
ment is very good (down to a few percent, compatible to the 
number of realizations). However, this agreement is not satisfac- 
tory for the non-diagonal elements when working under aggres- 
sive masks: in this case (and also when computing the dispersion 
of the statistic 0) we shall use the results obtained from 10,000 
Monte Carlo simulations. This assures a fair estimation of the 
correlation between different ap estimates and hence an accu- 
rate estimation of the overall statistic. 



3. Comparison to the Angular Cross Power 
Spectrum 

In this Section we shall compare the matched filter (as defined 
above) to the Angular Cross Power Spectrum (hereafter ACPS) 
method. This comparison will be made within the model moti- 
vated in the previous Section, 



t + am. 



(12) 



Note that, according to our notation, (a), = a, denotes the i-th com- 
ponent of the array a, where sip denotes the p-th array of some larger 
group of arrays. Same for matrices: (C),j denotes the array element in 
the !-th row and j-th column, not to be confused with Cp, which denotes 
the covariance matrix computed from arrays a,, and a,. 



and will be restricted to the large angular scales. This choice 
is motivated by the study of the Integrated Sachs-Wolfe (ISW) 
effect that follows in subsequent Sections of the paper, and that 
is typically restricted to / < 50. When studying small scales, one 
has to be careful with the SVD decomposition, which might fail 
for too large matrices. The matched filter in real space is, in most 
of those occasions, more adequate. 

The Angular Cross Power Spectrum (hereafter ACPS) can be 
viewed as a simplification of the matched filter presented here, 
where the covariance matrix is approximated by a diagonal ma- 
trix with identical non zero elements. Following the notation of 
Section im the estimate of a provided by this method is given 
by 



0!ACPS 



(13) 



where i/_,„ and m; ,„ are the Fourier multipoles of the signals s 
and m, respectively. Note that we shall refer to these signals in 
Fourier space, and thus the vectors s and m will contain the com- 
ponents of the 5/ ,„ and m/.,„ multipoles as explained in Section 



In order to compare this method to the matched filter, 
we have to define t and m and build s according to equa- 
tion (fT2l i. Throughout this paper, we shall not use real data 
but only Gaussian realizations generated from a given cosmo- 
logical model. For t, we choose CMB realizations for which 
the ISW contribution has been subtracted. I.e., we simulate 
the Fourier multipoles f/,m-s from an angular power spectrum 
computed using a modified version of the CMBFAST code 
(Seljak & Zaldarriaga(199^ with a cosmo logical parameter set 
equal to that given in Sper gel et al.(2007)[ = 0.1994, Oa = 

0. 759, Q.h = 0.0416, n, = 0.958, erg = 0.75 and t = 0.089. This 
will be the reference cosmological model hereafter. The template 
m is a Gaussian realization of a projection of the linear density 
field as computed from the matter power spectrum obtained with 
the same cosmological parameters. This density field is placed 
within a shell centered at z = 0.8 with a total width of Az ~ 0.8, 

1. e., the redshift range where ISW contribution is maximal (this 
will be addressed in detail in Section (14. It ). To each realiza- 
tion of t we added the component am (for a given choice of a, 
a - IQ^^), and applied the two methods. All maps were con- 
volved with a Gaussian PSF of FWHM= 2°. Let us remark that 
by this exercise we do not attempt to simulate ISW observations, 
but simply test the two methods in the context of equation ( fT2l i. 

In this comparison, we applied both the matched filter and 
the ACPS under two different masks: the first one assumes full 
sky coverage (fsky - 1), whereas the second one adopts the 
mask provided by the fourth data release of Sloan Digital Sky 
Survey (SDSS-DR4, Eisenst ein et al.(2001 ) ) combi ned with the 
KpO mask used in WMAP data analyses, ( [Hinshaw eFal.(2003)[ 
see left panel of Fig. ([T]i)- In Table ^ we display the results 
after applying both methods to an ensemble of 10,000 simula- 
tions. We consider two scenarios: (i) a unique /-bin limited to 
hnax = 15, / e [2, 15], and (ii) a set of 16 /-bins, ranging from 
/,„,„ = 2 to Imax - 50. We always find that the estimates of 
a are unbiased for both method^ and that the matched filter 
provides an estimate of the dispersion of a (see equation (|2]i) 
that actually agrees with the value recovered from the Monte 



' A mask in real space involves a convolution of different a values 
in Fourier space, which may generate a bias if a is not constant versus 
/.The case considered in this Section observes a constant a for every 
multipole, but in subsequent sections this will not be the case and a bias 
will appear. 
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Fig. 1. Three different masks used in this paper: Left: Mask corresponding to the product of the SDSS-DR4 mask times the KpO 
mask used in Hinshaw et al.(2003)| MiWte: Mask assigned to the future PAU-BAO galaxy survey times the KpO mask. Right: Mask 
of the NVSS times the KpO mask. 





hnax 1 ^ 






MF 1 ACPS 


MF 


ACPS 




(Q'>/tr,i 




(li)lo-i, 


(xl) 




All sky 


25.70 


16.50 


2,312 


192 


2,316 


192 


SDSS DR4 


34.68 


6.44 


271 


63 


201 


55 



Table 1. Under two different masks (SDSS-DR4 mask and all sky), we compare the performance of the matched fil- 
ter (MF) and the ACPS when trying to estimate the amplitude a from a data set s - t + am given s, m and 
the power spectrum of the Gaussian noise t, which is taken from a CMB power spectrum. In case (i) we con- 
sider a single /-bin containing all multipoles in I e [2,15], whereas in case (ii) we consider 16 Z-bins: / e 
[2, 3], [4, 5], [6, 8], [9, 14], [15, 25], [26, 28], [29, 31], [32, 34], [35, 37], [38,40], [41, 43], [44,45], [46, 47], [48, 49] and [50, 51]. We 
are quoting the results for the a statistic in case (ij, and for the;^'^ and y6 statistics in case (ii). Note that, in this case, the;^'^ statistic 
has been normalized by the number of degrees of freedom, i.e., the number of Z-bins. 



Carlo simulations. In case (i) we obtain that the matched fil- 
ter works better than the ACPS under the two masks consid- 
ered. Note that the noise signal in these analyses corresponds 
to the CMB (after having the ISW component subtracted), and 
that therefore the noise power spectrum scales as (|f/,,„P) oc l^^. 
On the other hand, the angular power spectrum of our density 
template scales roughly as oc const at / < 50. This 

means that for the low-Z range considered in ( i), the matched fil- 
ter is going to weight more the high-Z end multipoles: this effect 
makes this method superior to the ACPS (which weights all mul- 
tipoles equally) in the all sky case. For the SDSS-DR4 mask, the 
matched filter also accounts for the coupling among mutipoles, 
and this enlarges the difference between the two methods. One 
obvious question that arises is: how can the matched filter per- 
form better under the SDSS-DR4 mask than in the all sky case? 
The low I modes are degenerate under the SDSS-DR4 mask, that 
is, they are not orthonormal as for fsh - 1 and are decomposed 
onto other modes corresponding to smaller angular scales. I.e., 
the noisiest (low modes under the full sky mask are not eigen- 
modes anymore and they are partially dropped from the analysis 
(the matched filter method handles 176 different modes under 
the SDSS-DR4 mask, as opposed to 252 modes if f^ky - 1). 
This means that under the SDSS-DR4 mask we have a different 
statistic (since it handles a different number of degrees of free- 
dom) that is more concentrated in angular scales where the noise 
amplitude is smaller This provides this new statistic a better S/N 
ratio. 

In case (ii) and f^ky = 1, it turns out that, given the scaling of 
the power spectra of m and t, most of the S/N ratio is in the few 
higher I bins, centered at multipoles / - 45,46, ...,50. Indeed, 
these last multipoles are dominating the sums in equation ( fT3] l. 
For these few high-Z bins carrying most of the information, the 
change of the noise properties is very small, and therefore each 



of these bins is roughly equivalent to the rest. The weighting ap- 
plied by the matched filter introduces very little difference, and 
both methods perform similarly, yielding almost identical val- 
ues of x^- But again, under the SDSS-DR4 mask the coupling 
among multipoles is observed by the matched filter, and this in- 
troduces a difference between the two methods. However, it is 
clear that the matched filter proves comparatively better when 
when a single /-bin is considered, (case (i)). 

4. Application to ISW Studies 

The Integrated Sachs Wolfe (ISW) effect arises as a consequence 
of a late time variation of the gravitational potentials in the 
large scales. If there is a net change in the depth of the poten- 
tial wells while they are being crossed by CMB photons, then 
this radiation field will experience a gravitational red/blueshift. 
Crit tenden & Turok (1996)| pointed out that gravitational poten- 
tials should be traced by the Large Scale Structure (LSS), and 
proposed the cross-correlation of CMB maps to LSS surveys to 
unveil this signal. However, in most plausible models the time 
variation of the potentials occurs at late times (or low redshifts, 
z < 2), and the angle subtended by the linear scales for which the 
ISW effect is important is rather large (0 > 3 - 5°). This means 
that there will be room for relatively /ew independent ISW spots 
on the sky, i.e., this effect will be considerably limited by Cosmic 
Variance. Further, it is in this large angular range where the 
Galactic emission is more important, and errors in its subtraction 
might be more relevant. For this reason, it becomes necessary the 
use of masks that project out regions where this galactic contam- 
ination is large and cannot be removed accurately. Furthermore, 
when doing a cross-correlation analysis between CMB maps and 
LSS maps, the latter may not likely cover the whole sky, but also 
be restricted to a given limited region. In this context, it becomes 
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Fig. 2. (a) The dot-dashed line displays the cross power spec- 
trum of the projected density field through a shell centered at 
z - 0.8 (dashed line) times the ISW component (thin solid line). 
The thick solid line displays the total CMB angular power spec- 
trum, (b) Normalized signal to noise ratio (as given by equation 
fn} ) versus multipole /: the solid line corresponds to a density 
probe in a shell centered at z = 0.8, whereas for the dotted and 
dashed lines the central redshifts are 0.4 and 1.3, respectively. In 
both panels, the width of each density shell is equal to 20% of 
the comoving distance to the central shell redshift. 



particularly important the implementation of as sensitive as pos- 
sible cross-correlation tools that are able to handle optimally the 
limitations imposed by the sky masks. 

4.1. The S/N ratio in ISW Studies 

In this Subsection we briefly describe the ISW - LSS cross- 
correlation in the frame of the WMAP3 cosmogony. Unlike 



Douspis et al.(2008) we refrain from addressing the dependence 



of this correlation under different Dark Energy models. We con- 
centrate on a single LSS survey, and search for its optimal red- 
shift in terms of ISW detection. We study the amount of S/N 
that arises in those cases, and the angular scales where it is gen- 
erated. This sets the scenario for our cross-correlation method 
comparison. 

The ISW - LSS cross correlation arises from the fact that 
both LSS probes and gravitational potentials are tracing the un- 
derlying matter density field. The expression for the temperature 
anisotropies introduced by a gravitational blue/redshift reads 
( |Sachs & Wolfe(I967)|p'artinez-Gonzalez et al.(1990)| 



dt (j(fi, f). 



(14) 



The symbol (p denotes the time derivative of the gravitational 
potential (p. In linear theory, this expression can be rewritten as 
(e.g.,|Cooray(2002)|i 



JSW 



/ 



dr ji(kr) 



-^^mHl d(D/a) 
fe2 dr ' 



(15) 



The multipole coefficients a'^^^ are related to the ISW temper- 
ature anisotropies by equation ([3]). In this equation, r denotes 
comoving distance, k comoving wavevector, ji{x) denotes the 
spherical Bessel function of order /, Hq is the Hubble parame- 
ter, fl(r) is the scale factor and D(r) is the standard linear growth 
factor. The 3D Fourier mode of the density contrast is denoted 
by Note that for an Einstein-de Sitter Universe (D - a) the 
whole integral vanish. In a similar way, the multipole coefficients 
for the angular number density of a matter density probe (which 
will be taken to be galaxies in what follows) read 



af =(-0'(4: 



/ 



(2^ 



dr ji(kr) r^ ng(r)b(r, k) D(r) 6^, 



(16) 



with rtgir) the average number density of galaxies. The bias func- 
tion b(r, k) accounts for usual probes of the LSS actually being 
biased tracers of the underlying mass distribution (b > 1). This 
expression neglects the presence of shot (Poisson) noise in the 
galaxy number Note that the coordinate r here is being taking 
as a look-back time coordinate, equivalent to conformal time or 
redshift (z). For the sake of clarity, in what follows we shall use z 
as look-back time coordinate. We have that a given galaxy survey 
will probe the redshift range given by the product n(z) = b Ng, 
with Ng(z) = r^{z) ng(z)D{z). Note that, at this stage, we are ig- 
noring the k dependence of the bias function b. For simplicity, 
we shall rewrite equations (115116b as 



ISW,g 



= (-0'(4: 



with the A'^^'^(k, z) being referred to as transfer functions of the 
ISW and the galaxy fields, respectively. In real space, the cross 
correlation function between ISW temperature anisotropies and 
LSS probes reads 



dk 

(2^ 



X A 



ISW,g 



ik,z), 



(17) 



C 



Pi{cos0), 



with c'^^'^s the cross power spectrum. 



C 



IS W»g 

I 



k^dkA 



ISW 



Af P„,ik), 



(18) 



(19) 



and P,„(k) is the linear matter power spectrum. The symbol "®" 
denotes cross correlation. The theory makes actual predictions 
on this cross-power spectrum, and its covariance matrix is diag- 
onal if fsky - 1, and in general is simpler than that of the corre- 
lation function, (see the detailed analysis of |Cabre et al.(200"7)l l. 
As an example, we show a cross power spectrum (dot-dashed 
line) in the top panel of Fig. The total CMB contribution 
(for the chosen ACDM model) is given by the thick solid line, 
and its ISW component is displayed by the thin solid line. The 
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angular power spectrum corresponding to the projection of the 
galaxy field whose window function n(z) is centered at z = 0.8 
is displayed by the dashed line. The total width in redshift space 
is roughly Az ~ 0.80, so it is a thick shell. The cross power spec- 
trum peaks at scales at around / ~ 30 - 50, but its amplitude at 
/ ~ 200 is roughly equal to that at / ~ 2. This might suggest that 
the there is so much cross-correlation signal at large (/ < 10) as 
in small (/ > 100) scales, but indeed most (90%) of the signal 
comes from the large scales (/ < 40) if fskv = 1, as we will show 
below. Note that the amplitude of the density power spectrum 
(and hence the cross power spectrum) are taken arbitrary. The 
signal-to-noise (S/N) ratio for the measurement of a given mul- 
tipole / of the cross power spectrum can be computed once one 
takes into account that the aifl multipoles are real defined, and 
that the real and complex components of the ai^„ (m > 0) mul- 
tipoles, besides being equivalent and independent, must satisfy 
the constraint for the total amplitude (|fl/,mP) - Ci. We obtain 



[cf^^Cf + (Cf"'®«)2](Z/2+ 1)' 



(20) 



In this equation, Cf^^ is the CMB angular power spectrum and 
Cf is the LSS probe auto power spectrum. The quantity 



si = 



(21) 



is displayed versus I in the bottom panel of Fig. Q. I.e., this 
figure shows the ratio of the signal to noise ratio contained be- 
low some given I. The thick solid line corresponds to the galaxy 
survey centered at z = 0.8 mentioned above, whereas the dot- 
ted line for a galaxy survey with a window function centered at 
z = 0.4. The dashed line corresponds to a case where the galaxy 
survey is probed at z = 1.3. In all cases we are taking a shell 
width equal to 20% of the comoving distance to the peak of the 
window function 11, and we are assuming that fs^y - 1 . We see 
that regardless where Tl(z) peaks, practically half of the total sig- 
nal is contained at multipoles I < 10, whereas its 90% fraction 
is typically contained at I < 40 - 50. (Had we considered thin- 
ner shells -width equal to 2% of the comoving distance-, then 
all those shells below z = 0.8 would have still shown a pat- 
tern very close to that given by the solid line). This suggests 
that by dropping all multipoles above / - 50 (or by neglecting 
scales smaller than 6 « 3°- 4°) one should recover practically 
the same ISW detection significance. This sets a useful con- 
sistency check, given the number of other physical phenomena 



(Rees-Sciama effect ( Rees & Sciama(1968)), kinetic Sunyaev- 
Zel'dovich effect ( jSunyaev & Zeldovich(1972)) , intrinsic source 
emission, etc) that arise at smaller angular scales and that, a pri- 
ori, correlate with the spatial position of LSS probes. 

In Fig. (O we show how the total signal to noise ratio de- 
pends on the central redshift for the galaxy survey window func- 
tion n(z). For the thick solid line, the width is taken to be roughly 
20% of the comoving distance to the central redshift. In redshift 
space, it implies a width of Az ~ 0.7 for central redshift z = 0.3, 
Az ~ 0.84 for central redshift z = 0.8, and Az ~ 0.99 for central 
redshift z = 1.3. The thin solid line observes a width of only 2% 
the comoving distance to the central redshift, and this translates 
into Az ~ 0.07, 0.08 and 0.1 for central redshifts z = 0.3, 0.8 
and 1.3 respectively. We should obtain the larger detection lev- 
els for thick shells (large Az's) and z - 0.8, and this motivates 




1.0 1.5 2.0 
Redshift z 



Fig. 3. Total signal to noise ratio of the cross power spectrum 
detection as given by the numerator of equation ( |2TI ) with re- 
spect to the central redshift of the density probe shell. The thick 
solid line displays the case when the shell containing the density 
probes has a width equal to 20% of the comoving distance to the 
shell. For the thin solid line, this width is only 2%. 

our choice of a thick survey centered at this redshifQ. Note that 
our assumption that the bias is independent of scale might not be 
accurate, but it has less impact in the large scales (low Is) where 
most of the effect is coming from. We do not expect significant 
changes after introducing a scale dependent bias in our galaxy 
survey description, although we shall address this issue in detail 
when applying our method to real CMB and LSS data. Note that, 
a priori, this method can be applied the same on multiple redshift 
shells, and is affected in exactly the same way than the ACPS by 
realistic aspects such as the redshift or scale dependence of the 
bias, survey incompleteness, etc. 

We next compare the performance of the matched filter to 
that of the ACPS. We use one Gaussian realization of our chosen 
density shell, and compute a single Gaussian realization of an 
ISW map compatible to it. If a^'^^^ are the Fourier multipoles of 
our density 2D template, then they can be related to those of a 
compatible ISW map via (e.g., |Cabre et al.(2007)[ l 



JSW 



aia^ 



c 



l,m 



(22) 



The Gaussian signal /?/ ,„ is the part of the ISW component that 
is uncorrelated to the LSS 2D template, verifying (|/S/,,„P) = 
C'^^ -(C'l^^'^^)^ /C^ , where denotes the angular power spec- 
trum of the LSS probe. Note that this is a correct way to ex- 
press the ISW field in terms of the galaxy density field as long 
as both fields are Gaussian and their are completely determined 
by the first and second order momenta. The ratio c['^*^*^/C* is 
explicitly identified with ai, which is precisely the output of our 
matched filter technique. According to the theory, c; shows a 
strong dependence on /, and therefore our method must be ap- 
plied in separate /-bins. Gaussian simulations of the CMB were 
built from the addition of om fixed ISW template plus realiza- 
tions of a CMB angular power spectrum for which the ISW 
component had been subtracted, just as for the t component 



* It has been noted elsewhere (e.g. Afshordi(2004) I that by combin- 
ing different LSS surveys at different redshifts one can obtain larger 
S/N ratios. We shall avoid that discussion here and focus our method 
comparison on one single survey. 
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in Section (O. The realizations from the modified CMB angu- 
lar power spectrum were computed upto a maximum multipole 
I - 160, and convolved with a Gaussian beam of 2°of FWHM. 
The fixed ISW and the LSS maps were also convolved with the 
same PSF, and all maps were produced under the HEALPijfl 
( |G6rski et al.(2005)l ) resolution parameter Nside = 64. 

4.2. Performance under Different Masks 

The simulated maps were cross-correlated to the projected 
density map (hereafter denoted as m) with both the ACPS 
and the matched filter methods, according to the multipole 
decomposition given in Section 12.21 Their performance was 
compared under three diff'erent masks shown in Fig. ([TJ: the 
left panel shows the mask corresponding to the sky cover- 
age of the fourth data release of the Sloan Digital Sky Survey 
(SDSS-DR4, Eis enstein et al.(2001)| l. This mask was multi- 
plied by the KpO mask used in the analysis of WMAP data 
( Hinshaw et al.(2003)| l, and therefore the combined mask ob- 
serves a bit less than 10% of the total sky. In the middle panel we 
consider the fraction of the sky covered by the upcoming survey 
PAU-BAO E This survey is planned to cover ~ 10,000 square 
degrees of the celestial northern hemisphere, and in this work 
we have assumed that it is limited to the region b > 20° outside 
the KpO mask, in such a way that fsi^y - 0.26. Finally, the right 
hand side panel displays the product of the KpO mask with the 
mask corresponding to the NVSS survey ( [Condon et al.(199"8)l l. 
In this case, fsky - 0.65. 

A total of 10,000 simulations were run in order to perform 
the method comparison, and results are given in Table (|2]). The 
sensitivity of both methods is measured by the ™d the yS 
statistics for different choices of the maximum multipole l^ax 
considered in the analyses. We see that, unlike in the previous 
Section (where the matched filter was in general significantly 
more sensitive under the data model s = t H- am), in this ISW 
context both filters perform very similarly. There seems to be 
however a slighter better sensitivity of the matched filter under 
the most aggressive masks, but the difference is small (at least 
in terms of the output of the;^'^ and statistics). We remark that 
the ACPS is the Legendre transform of the angular cross cor- 
relation function, and that both methods should, a priori, show 
similar sensitivities. Of the two statistics considered in Table (|2|i, 
y6 provides the largest significance of the ISW detection (its dis- 
tribution is very close to Gaussian, and the number of sigmas 
yield smaller chance probabilities), but is remarkable that for 
both of them this significance finds a maximum at Z ~ 5 -40, and 
then starts dropping again. According to this result, by merely 
observing multipoles below I — 40 we should recover practi- 
cally all the ISW-LSS cross correlation significance. This points 
in the direction of our si estimates in Subsection ( 14.1b . which 
predicted most of the S/N ratio to be confined at low Z-s, and 
should prevent confusion with other secondary effects (giving 
rise to correlations of other nature such as Rees-Sciama effect, 
kinetic Sunyaev-Zel'dovich effect, etc). 

One advantage of the correlation methods implemented here 
is that they conduct analyses in separate ranges of angular 
scales, enabling a direct comparison with theoretical predictions. 
This is explicitly shown in Figs. (|4]i and (|5]l. In all panels of Fig. 
(HI, solid lines display the theoretical prediction for the correla- 
tion coefficient a; = c'^^'^^ /C^ versus multipole /. Our choice 
for the density and ISW templates is such that full sky analyses 



should yield a/ estimates exactly on the solid lines. Filled circles 
and triangles display average estimates of or; for the matched 
filter and the ACPS methods, respectively. Error bars denote 
the rms scatter for each of them. Let us first note that, accord- 
ing to the error bars displayed in Fig. (|4]l , most of the infor- 
mation seems again to be restricted to the large angular scales 
(/ < 30 - 40), as it has been quoted above. 

Let us now address the issue of the impact of the mask on 
the methods' output. A clear low bias can be seen in the es- 
timates of Of; for low l-s, specially under the SDSS-DR4 and 
PAU-BAO masks. This is a direct effect of the mask: a multi- 
plication of the full sky map by the actual mask in real space 
translates into a convolution in Fourier space, which involves a 
wider range of multipoles the smaller the mask is. Therefore, 
the ai estimates for the SDSS-DR4 mask will be the result of an 
average or smoothing of the a/ values in a wide space of I mul- 
tipoles. Since at low l-s, the values of o-/ are falling steeply, the 
convolution will provide a value that is smaller than the actual 
theoretical value, as it is displayed in the left panel of Fig. (HJ. 

At the same time, the mask introduces another bias in the 
large I range, for which the actual recovered values of ai are 
above the theoretical expectation. This is showing the fact that 
the aliasing introduced by the mask is shifting some large scale 
(low /) power into the small (large I) angular range. We at- 
tempted to quantify this aliasing by performing the following 
exercise: we generated one ISW map by using multipoles re- 
stricted to the range I e [2, 10]. We multiplied this map by the 
SDSS-DR4 mask used in this work, and computed the power 
spectrum of the resulting map. We measured the aliased vari- 
ance contained above some Z,„,„ > 10 by using cr^ [/„,,„] = 
Y,i„J2l + l)/(47r)C/. We found that about 30% of the total rmfl 
was contained above = 30. I.e., the ISW power aliasing from 
large to small scales is indeed significant. Let us remark as well 
that this effect is more present for the matched filter ai estimates, 
as we shall discuss next. 

A direct visual comparison of the two methods can be found 
in Fig. (Is), where the S/N ratio for each multipole bin is shown. 
The matched filter (solid circles) performs more accurately than 
the ACPS method (filled triangles), specially under the SDSS- 
DR4 and PAU-BAO masks, although this has a limited impact on 
the final detection significance quoted by the and statistics 
(see next Section). 



5. Discussion and Conclusions 

In order to assess the sensitivity of the two methods to the ISW, 
we have defined two different statistics: the statistic uses a 
quadratic combination of the method's output (in a similar way 
Tegmark(1997)[ l, and the yS statistic, which instead is Un- 



as m 

ear in the a/'s and for our purposes is Gaussian distributed. Both 
statistics pick up the information of the cross-correlation in dif- 
ferent ways. The;^'- statistic is more sensitive to the presence of 
ISW at the very large angular scales only, and rapidly gets de- 
graded as smaller angular scales are considered, i.e., it seems to 
be particularly affected by the inclusion of modes that have a low 
S/N ratio. On the other hand, the statistic is more sensitive to 
the actual signal to noise ratio even at scales where such ratio is 
below unity, and, as mentioned above, seems to be more efficient 
in terms of detection of the ISW-LSS cross-correlation. 

This connects to the multipole or range where most of the 
correlation is arising. Of the two methods, the matched filter 



http://www.healpix.jpl.nasa.govJ ' The total rms was computed by taking /„,;„ = 2, and was less than 

http : //www . ice . csic . es/research/PAU/PAU-welcome .html | 10% off the estimate obtained from the map in real space. 
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/ - 31 


hnax 40 








MF 


ACPS 


MF 


ACPS 


MF 


ACPS 


MF 


ACPS 


MF 


ACPS 


SDSS-DR4 


$)l<rB 


1.79 
1.24 


1.31 
0.65 


2.04 
1.87 


1.56 
1.35 


1.61 
1.74 


1.48 
1.70 


1.46 
1.57 


1.40 
1.61 


1.33 
1.46 


1.30 
1.22 


PAU-BAO 




2.93 
1.69 


2.02 
1.28 


2.94 
2.48 


2.45 
2.22 


2.12 
1.93 


2.29 
2.33 


1.87 
1.90 


2.02 
2.33 


1.61 
1.69 


1.74 
1.98 


NVSS 




6.27 
3.24 


7.03 
3.37 


5.30 
3.51 


6.04 
3.70 


4.15 
3.44 


4.71 
3.62 


3.47 
3.43 


3.91 
3.56 


2.75 
3.19 


3.06 
3.19 



Table 2. Comparison of the matched filter (MF) to the ACPS in the context of ISW studies. We quan- 
tify the sensitivity of each method by two statistics: (which has been normalized by the num- 
ber of degrees of freedom) and p, for different choices of Imax- In total, we considered 21 Z-bins: I e 
[2,3],[4,5],[6,8]J9a4],[15,25],[26,28],[29,31],[32,34]J35,37],[38,40],[41,43],[44,45],[46,47],[48,49],[50,51],[52,53],[54^^ 
and [60,61]. 




seems to be more confined in /-space than the ACPS when look- 
ing at the output of the p statistic: it quotes the maximum signifi- 
cance at l„ax - 14 and always drops at larger l,„ax-^, whereas the 
ACPS seems to peak at around l„jax = 30 - 40. In this case, the 
exception is the NVSS-like survey, for which fi yields the max- 
imum detection significance at l,„ax - 14. This different behav- 
ior suggests that aliasing induced by SDSS-DR4 and PAU-BAO 
masks is indeed shifting some S/N into the smaller scales (larger 
/-s), but this effect is not perceptible beyond l,„ax = 30 - 40. This 
is also visible in Figs.(|4]and|5]i: under the most aggressive masks, 
there is less information in the first Z-bins {I e [2, 3], [4, 5]), 
whereas for the NVSS-like survey these contain the largest val- 
ues of the S/N ratio. 

These two figures also show that, in almost every /-bin, the 
matched filter ai estimates seem to be more accurate than the 
ACPS method, and that, at the same time, they seem to be more 
affected by aliasing. These two facts are connected: the matched 
filter tends to pick up the signal from those modes having largest 
S/N ratio within each Z-bin. For small fsky and moderately high 
/-bins, these modes are actually aliased components of low I 
modes whose power has been shifted by the mask into smaller 
scales. This makes the matched filter provide more accurate a/ 
values in these high /-bins, but these estimates are actually highly 
correlated to those found at lower /-s. This limits the amount of 
information that high-/ bins actually add, and partially explains 



the high bias of the a/ estimates at large /-s in the left panel of 
Fig.®. 

In this work, we have generalized the implementation of the 
matched filter into the Fourier space of the 2D sphere, and ap- 
plied it in the context of CMB analyses and ISW studies. The 
matched filter provides a tool to estimate the level of presence 
of some template m in some measured signal s containing a 
noise component t. This tool is optimized for the case of t be- 
ing isotropic and Gaussian distributed, and hence is particularly 
suited for cross-correlation tests where the CMB is the back- 
ground (noise) signal. In Fourier (or multipole) space, the corre- 
lation properties of the CMB are particularly simple (specially, 
but not only, in the full sky case fskv - 1). For fsky < 1, the 
covariance matrix of the CMB multipoles can be inverted via 
a SVD approach: this permits simultaneously identifying those 
Fourier modes containing more information and dropping those 
other modes that introduce numerical error. After all modes have 
been sorted in terms of their S/N ratio, the matched filter al- 
gorithm weights them accordingly in order to produce an op- 
timal (minimum variance) output for the cross-correlation test. 
We have compared this method with the standard Angular Cross 
Band Power Spectrum (ACPS), and found the the matched filter 
to be either superior or equivalent to the ACPS. 

In the context of ISW analyses, the matched filter provides 
estimates of the level of cross-correlation of CMB maps with 
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Fig. 5. Recovered signal to noise ratio for each multipole bin with the matched filter (filled circles) and the ACPS (filled triangles), 
under the three masks considered in Fig. (|4]i. 



LSS probes at separate multipole ranges, and this enables a di- 
rect and clean comparison to theoretical predictions . When ap- 
plying both the matched filter and the ACPS methods to three 
mock surveys having distinct values of fsky, we find that both 
methods perform similarly (the matched filter is slightly more 
sensitive under aggressive masks, the ACPS more accurate un- 
der the NVSS mask). The masks introduce some power aliasing 
from large into small angular scales, but this does not prevent 
most of the S/N ratio of the ISW-LSS cross correlation from 
being confined into the large angular scales (I < 40). This l- 
confinement may result particularly useful when distinguishing 
this effect from other secondary anisotropies that, while tracing 
the LSS distribution, arise at smaller angular scales. 

Natural extensions of this work involve large angle compo- 
nent separation in future CMB maps when tracers or templates 
for the signal to be distinguished are available, such as galactic 
or extragalactic radio, synchrotron or dust maps, large scale 1/f 
noise component, local kSZ or tSZ contributions, etc. 

Acknowledgements. I am grateful to Tony Banday for carefully reading the 
manuscript. I acknowl edg e the use of the HEALPix jGorski et al.(2005)t pack- 
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